fitL_explicit = glm(y~age, family=gaussian(link="identity"), data=d)

fitL <- glm(accuracy ~ age, data = d)

effL <- data.frame(
  allEffects(
    fitL,
    xlevels = list(age = seq(min(d$age), max(d$age), .05))
  )[["age"]]
)

ggplot(d, aes(x = age, y = accuracy)) +
  coord_cartesian(ylim = c(0, 1)) +
  geom_point(size = 4, alpha = .5, color = "darkblue") +
  geom_ribbon(
    data = effL,
    aes(x = age, ymin = lower, ymax = upper),
    alpha = .3, fill = "darkred", color = NA,
    inherit.aes = FALSE
  ) +
  geom_line(
    data = effL,
    aes(x = age, y = fit),
    size = 2, color = "darkred",
    inherit.aes = FALSE
  ) +
  theme(text = element_text(size = ts, color = "black")) +
  scale_x_continuous(breaks = seq(floor(min(d$age)), ceiling(max(d$age)), .5)) +
  scale_y_continuous(breaks = seq(0, 1, .1)) +
  ylab("accuracy") + xlab("Age (years)")

effLogit <- data.frame(
  allEffects(
    fitLogit,
    xlevels = list(age = seq(min(d$age), max(d$age), .05))
  )[["age"]]
)

p1 <- ggplot(d, aes(x = age, y = accuracy)) +
  coord_cartesian(ylim = c(0, 1)) +
  geom_point(size = 4, alpha = .5, color = "darkblue") +
  geom_ribbon(
    data = effLogit,
    aes(x = age, ymin = lower, ymax = upper),
    alpha = .25, fill = "blue", color = NA,
    inherit.aes = FALSE
  ) +
  geom_line(
    data = effLogit,
    aes(x = age, y = fit),
    linewidth = 2, color = "blue",
    inherit.aes = FALSE
  ) +
  theme(text = element_text(size = ts, color = "black")) +
  scale_x_continuous(breaks = seq(floor(min(d$age)), ceiling(max(d$age)), .5)) +
  scale_y_continuous(breaks = seq(0, 1, .1)) +
  labs(y = "accuracy", x = "Age (years)")

p1

d$group <- as.factor(d$group)

eff_obj <- allEffects(
  fitLint,
  xlevels = list(age = seq(min(d$age), max(d$age), .05))
)[["age:group"]]

effL <- as.data.frame(eff_obj)

if (!("fit" %in% names(effL))) {
  pred_col <- intersect(c("fit", "response", "fitted", "yhat"), names(effL))
  if (length(pred_col) == 0) stop("Nessuna colonna di predizione trovata in effL.")
  effL$fit <- effL[[pred_col[1]]]
}

ggplot(d, aes(x = age, y = accuracy, shape = group, color = group)) +
  geom_point(size = 4, alpha = .6) +
  geom_ribbon(
    data = effL,
    aes(x = age, ymin = lower, ymax = upper, fill = group, group = group),
    alpha = .3,
    color = NA,
    inherit.aes = FALSE
  ) +
  geom_line(
    data = effL,
    aes(x = age, y = fit, color = group, group = group),
    linewidth = 2,
    inherit.aes = FALSE
  ) +
  scale_color_manual(values = c("darkorange3", "darkgreen")) +
  scale_fill_manual(values = c("darkorange3", "darkgreen")) +
  theme(text = element_text(size = ts, color = tc)) +
  scale_x_continuous(breaks = seq(floor(min(d$age)), ceiling(max(d$age)), .5)) +
  scale_y_continuous(breaks = seq(0, 1, .1), limits = c(0, 1)) +
  labs(y = "Accuracy", x = "Age (years)")

fitLogitint <- glm(
  accuracy ~ age * group,
  data = d,
  family = binomial(link = "logit"),
  weights = rep(k, nrow(d))
)

effP <- data.frame(
  allEffects(
    fitLogitint,
    xlevels = list(age = seq(min(d$age), max(d$age), .05))
  )[["age:group"]]
)

effP$group <- as.factor(effP$group)

ggplot(d, aes(x = age, y = accuracy, shape = group, color = group)) +
  geom_point(size = 4, alpha = .6) +
  geom_ribbon(
    data = effP,
    aes(x = age, ymin = lower, ymax = upper, group = group, fill = group),
    alpha = .3, color = NA,
    inherit.aes = FALSE
  ) +
  geom_line(
    data = effP,
    aes(x = age, y = fit, group = group, linetype = group, color = group),
    linewidth = 2,
    inherit.aes = FALSE
  ) +
  scale_color_manual(values = c("darkorange2", "darkgreen")) +
  scale_fill_manual(values = c("darkorange2", "darkgreen")) +
  theme(text = element_text(size = ts, color = "black")) +
  scale_x_continuous(breaks = seq(floor(min(d$age)), ceiling(max(d$age)), .5)) +
  scale_y_continuous(breaks = seq(0, 1, .1), limits = c(0, 1)) +
  ylab("accuracy") + xlab("Age (years)")